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ABSTRACT 



Context. The Sun's polar fields and open flux around the time of activity minima have been considered to be strongly 

correlated with the strength of the subsequent maximum of solar activity. 

Aims. We aim to investigate the behavior of a Babcock-Leighton dynamo with a source poloidal term that is based on 

the observed sunspot areas and tilts. In particular, we investigate whether the toroidal fields at the base of convection 

zone from the model are correlated with the observed solar cycle activity maxima. 

Methods. We used a flux transport dynamo model that includes convective pumping and a poloidal source term based 

on the historical record of sunspot group areas, locations, and tilt angles to simulate solar cycles 15 to 21. 

Results. We find that the polar fields near minima and the toroidal flux at the base of the convection zone are both 

highly correlated with the subsequent maxima of solar activity levels (r = 0.85 and r — 0.93, respectively). 

Conclusions. The Babcock-Leighton dynamo is consistent with the observationally inferred correlations. 

Key words. Magnetohydrodynamics (MHD) - Sun: dynamo - Sun: surface magnetism 



1. Introduction 

The solar magnetic cycle is believed to be the result 
of a magnetohydrodynamic dynamo. One family of these 
dynamo models is the Babcock- Leigh t on (B L) model, 
which was first propose d by Babcoclg (|1961h and fur- 
ther elaborated on by iLeightonl (|l964f ). The essence 
of the BL mechanism is that the generation of the 
poloidal field comes from the decay of the tilt sunspot 
groups on the solar surface. BL models have been suc- 
cessful in repro ducing some cha r acteri st ics of the so- 
lar cycle (e.g., Choudhuri et al. _ 19951: Durnevl 1997 ; 



Choudhuri fc Dikpatilll999t iDikpati fc Charbonneaulll999t 



Nandy fc Choudhurill2002HChatteriee et alJJ2004 iRempe 
2006t IGuerrero fc de Gouveia D al Pino 2007; Ye ates et al 
_ , including the solar cycle irregularities (e.g., 
Charbonneau fc Dikpatil 20001: Durnev [2000: Cha rbonneaul 



m 

ii). 



2001; Karak| l2010t IKarak fc Choudhurill201lT ) . Recent anal- 
yses of the long-term sunspot observations strengthen the 
idea that the dynamo is c o nsistent with th e BL m echanism 
(iDasi-Espuig et all. 1201(1 ICameron et all . 120101 [hereafter 
CJSS10]: lKitchatinov fc OlemskovL l20ilT 



The important ingredients in the BL dynamo model are 
(i) the generation of poloidal flux due to the emergence of 
tilted sunspot groups and the subsequent evolution of the 
surface magnetic flux, (ii) the transport of poloidal flux 
to the tachocline, (in) the generation of toroidal magnetic 
flux due to the winding up of the poloidal field (the $7- 
effect), and (iv) the subsequent formation and rise of flux 
loops through the convection zone until they emerge at the 
surface. This last process also includes the tilting of the flux 



tubes with respect to the equator. For more info r matio n 
about the models, see the review by ICharbonneaul <|2010T> . 

Some of these processes are directly observable, in par- 
ticular, the surface evolution of the emerged magnetic 
field, part (i). Empirically the large-scale evolution is de- 
scrib ed well by the surface flux transport (SFT) model 
(e.g..lWang et allll989|: ISchriivenl2PM IMackav et al.ir2 002: 



i Baumann et al.l 12004 Ijiang et alj 120101: iMackav fc Yeatesl 
I2012H . We previously (CJSS10) used the SFT model to 
reconstruct the surface field and open flux for the pe- 
riod 1913-1986. The observed sunspot longitudes, latitudes, 
areas, and cycle-averaged tilt angles were used to cre- 
ate the source term. The results of that model compare 
well with the op en flux derived from geomagnetic indices 
(Lockwood 2003) and with the reversal times of the polar 
fields (|Makarov et al.ll2003l) . 



The remaining processes (ii), (Hi), and (iv) deal with 
the subsurface dynamics, which are currently poorly con- 
strained by observations. These processes are not in- 
cluded in the SFT model but are described by the 
flux transport dynamo (FTD). The surface fields result- 
ing from the two models are consistent, provided the 
FTP model has sufficient near-surface convective pump- 
ing (jCameron et al.| |2012f). The e ff ect o f pumping has 
been studied by IKarak fc Nandvl (|2012h . who demon- 
strate that the effect of the convective pumping re- 
duces the memory of the solar activity to one cy- 
cle. Other effects of the pumping on the BL d ynam o 
have also been studied by, e.g . . iKapvla et al.l (I2 006). 
IGuerrero fc de Gouveia Dal Pinol (12 008). D o Cao fc Bmnl 
(|201ll ). and lKitchatinov fc Olemskovl (|2012t ). 
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Our approach is to study the Babcock-Leighton dy- 
namo model using as many observational constraints as 
possible. We thus take the historical RGO/SOON sunspot 
record as the basis for constructing the source of poloidal 
flux, we employ boundary conditions where the magnetic 
field is vertical at photosphere, and we impose substantial 
downwards pumping near the s urface (as required by ob- 
servations, Came r on et al.|[2012f ). The differential rotation 
(jSchou et al.lll998|) , near-surface meridional flow, and near- 
surfa ce turbulent diffusivity (see references in ICameronl 
|2011|) are observationally constrained and modeled in a sim- 
ilar way. 

The diffusivity and pumping throughout the bulk of 
the convection zone, as well as the deeper structure of 
the meridional circulation remain observationally uncon- 
strained. In this paper we first construct a reference case 
with intermediate choices (between the extremes of what 
can be considered as reasonable) for both the pumping and 
the diffusivity in the bulk of the convection zone. After pre- 
senting the results of the calculations for this reference case, 
we then consider the effect of varying these parameters. A 
more extensive study of the effect of the diffusivity and the 
meridiona l circulation in transporting the field is given by 
lYeates et al.l (J2008T) for the case without magnetic pump- 
ing. The meridional circulation deep below the surface is 
controversial, with opposing views on such basic questions 
as whether there are one or two cells in the radial direction. 
For this paper we have assumed a one-cell structure. 

Here we concentrate on processes (i), (ii), and (Hi), us- 
ing the FTD model to study the evolution of the toroidal 
flux when the source of poloidal flux, process (i), is based 
on the historical sunspot record. In a full BL dynamo model 
process (i), the appearance on the surface of tilted sunspot 
groups, is intimately connected with the formation and 
rise of flux loops from the bottom of the convection zone, 
i.e., process (iv). Since we are using the observed sunspot 
record for process (i), we can omit the formation and rise 
processes entirely. Ignoring processes (iv) means that we 
are solving a driven system rather than a self consistently 
driven dynamo. However, it allows us to address the ques- 
tion of whether the polar field strength at the minimum 
produced by a BL dynamo with a poloidal source term de- 
rived by observations of sunspot groups and cycle-average 
tilt angles is correlated with the observed strength of the 
foll owing cycle. To be co nsistent with observations (e.g. 
see I Wang fc Sheelevll2009l ). a strong correlation should be 
present in the model. 

The correlation between the minima of the open flux 
(which is closely related to the Sun's axial magnetic dipole 
moment) and th e strength of the subsequent maxima o f 
solar activity Csee lWang fc Sheelevll2009HJiang et al.ll201ll ) 
implies processes (ii) and (Hi) are - on average over a cycle- 
essentially linear. The argument is that a nonlinearity in go- 
ing between the polar fields and the toroidal fields at the 
base of the convection zone would manifest itself in the 
number of sunspot groups that appear on the surface. Such 
a nonlinearity would then result in maxima of activity lev- 
els that are not strongly correlated with the polar fields. 
The maxima of the toroidal field at the base of the convec- 
tion zone, which is determined by processes (i-iii), should 
therefore also be strongly correlated with the maxima of 
the activity level. This is an observationally inferred result, 
and it is therefore a test for the model. 



There have been two previous efforts to relate the 
amount of toroidal flux produced by an FTD model and 
the observed suns p ot re cord. The first effort, started by 
iDikpati fc Gilmanl (J2006I ). used the time evolution of the 
area coverage of sunspots (a measure of the amplitude of 
the solar cycle) as input. A very idealized model was used to 
convert this into the source of poloidal flux. In particular, 
their source term for cycle n requires knowing of the timing 
of the minimum between cycles n — 1 and n. Their finding 
was that the toroidal field at the base of the convection zone 
during any cycle n in the model was highly correlated with 
the actually observed strength of cy cle n. An explanation 
for th is correlation is suggested by ICameron fc Schiisslerl 
(|2007t ) : incorporation of the timing of the minimum in the 
source term for th e poloidal flux, in conjuncti on with the 
Waldmeier effect (jCameron fc Schiisslerl 120081) introduces 
similar correlations even when the amplitudes are indepen- 
dent random realizations. 



The second effort is that of I Jiang et all (|2007l ). They 
found in the high-diffusivity regime that the poloidal flux 
is transported to the base of the convection zone at the 
same time as it is transported to the poles. The simulta- 
neous transport of the field to the poles and to the base of 
the convection zone presented in this model offers an expla- 
nation for the observed correlation between the polar fields 
and the toroidal field at the base of the convect ion zone (see 
also lChoudhuri et al.ll2007t ). lJiang et al.l (J2007h then scaled 
the poloidal field at the end of each cycle so that the polar 
fields of the model had the same strength as those observed 
by the Wilcox Solar Observatory. The winding up of this 
poloidal field by differential rotation produces a toroidal 
field that is correlated with the polar fields. The high cor- 
relation between the toroidal fields of the model and the 
observed level of activity then follows from the previously 
reported correlation between th e observed po lar fields and 
the amplitude of the next cycle (jSchatten et al. 1978). The 
exact nature of the sources in this model is irrelevant, since 
the field is rescaled, and this model says little about the 
causes of the cyclic amplitude variations. 

These previous attempts were carried out before 
cycle-to-cycle tilt angle variations had been reported by 
iDasi-Espuig et all (J2010I ) and before it was recognized that 
a vertical field boundary condition at the photosphere and 
strong pumping are required to make the FTD model match 
the observations of th e surface evolution of magnetic flux 
(jCameron et al.l l2012h . The question thus remains open 
as to what happens when we use the FTD model with 
these new constraints. H ere we address this q uestion using 
the FTD with pumping (jCameron et al.ll2012h and sources 
based on as much information as is available (see CJSS10) 
to follow the evolution of the subsurface magnetic field over 
cycles 15 to 21. 

The paper is organized as follows. The BL dynamo 
model is described in Section 2. The results for a reference- 
case set of parameters are presented in detail in Section [3l 
The influence of the observationally poorly constrained pa- 
rameters describing the magnetic diffusivity and pumping 
in the bulk of the convection zone is studied in Section |4] 
The conclusions and discussion of our results are given in 
Section [5) 
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2. Flux transport dynamo model 

The FTD model is based on the induction equation for an 
azimuthally symmetric field: 



8 A 1 1 



dB 
~dt 



where 



V {V 2 



• )B 



dr 



1 dr) d 
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(2) 



s(B p ■ V)ft, 



B = B(r,9)e^ + V x [A(r,9)e^] 



(3) 



Here B(r, 9) is the toroidal component of the magnetic field 
and A(r, 9) the toroidal component of the vector potential, 
related to the poloidal component of the magnetic field by 
B p = V x (Ae,p); v' is the sum of the meridional velocity 
field v(r,9) and magnetic pumping -y = 7(r)e r ; 0(r, 9) is 
the differential rotation profile; S(r, 9, t) is the source term 
for the poloidal field; r\ (r) is the magnetic diffusivity; and 
s = r sin (9. 

The code used to solve the above problem was developed 
at the MP S. It has be e n che cked against the benchmark 
dynamo of Jouve et all ([2008D and was previously used in 
ICameron et al.l ( 20121) . In essence it treats the advective 
term explicitly and an alternating direction implicit scheme 
for the diffusive terms. We use 181x71 grid cells in the lat- 
itudinal and radial directions, respectively, and a timestep 
of one day. 

2. 1 . Boundary conditions 

We carry out our calculations in a spherical shell (0.65-R© < 
r < Rq, < 9 < it). At the poles (9 — 0, it), we have 

A = 0, B = 0. (4) 

The inner boundary matches a perfect conductor, 

d{rB) 



A = 0, 



dr 



= at r = 0.65i? 



©• 



(5) 



As described in ICameron et al.l (J2012), the appropriate 
outer boundary is that the field is vertical there, 



d_ 

dr 



rA) = 0, B = at r = R e 



(G) 



Th e vertical outer boundary c onditi on is also proposed 
by Ivan Balleeooiien fe Mackavl (|2007l ) and has previously 
been used in the BL dynam o models by, for exam ple 



We also follow [Cameron et al.l (|2012l ) for the turbulent 
diffusivity profile 



V(r) = Vrz + 



Vcz ijrz 

2 

Vs - Vcz 



1 + crf 
1 + erf 



r — r r 
d~ 
r — r. 



(7) 



where the subscript rz is for the radiative zone, cz is for 
the convection zone, and s is for the surface properties, 
rj rz = 0.1 km 2 s" 1 , rj cz — 10 km 2 s _1 , ij s = 250 km 2 s _1 , 
r cz = O.7i?0, r s = O.95i?0, and d — 0.02i? Q . The argument 
for the choice of the surface diffusivity t] s is given in CJSS10 
- it is near the middle of the range give n by observations 
(see Table 6.2 of lSchriiver fc Zwaanll2000l) . This value for t] s 
is also consistent with the range found using com prehensive 
photospheric simulations (jCameron et al.l [20111 ). The tur- 
bulent diffusivity in the convective zone r\ cz is not directly 
constrained. The effect of different values for this parameter 
will b e studied in Section B[ 

In lCameron et all (|2012h . constraints were placed on the 
amount of magnetic pumping near the top of the convection 
zone. However, the strength of the pumping in the deeper 
layers is unconstrained. We assume here that the magnetic 
pumping extends to the base of the convection zone, albeit 
at a weaker level than in the near-surface layer. The radial 
dependence of the pumping is assumed to be 



7. 



7(rj = — 



where "f c 



1 + erf 

7s -7 



1 + erf 



(8) 



-2 m s _1 , 



Is 



-20 m s" 



0.9Rr. 



©) 



Deen used in tne hl oynamo model s by, tor exam ple, (jHoward et al.N1984t IHowa rd 1991; Sivaraman e t all 
iMunoz-Jaramillo et al.l (|2009| ) and iNandv et all (|2011l) . are averaged over each cycle as in lDasi-Espuig et alf ( 



and r 1 — 0.7Rq. We assume here that the enhanced near 
surface pumping extends slightly deeper than the region of 
enhanced surface turbulent diffusivity (r ns — 0.9Rq rather 
than r s — 0.95Rq), below which its strength drops to 
2 m/s. Since the value of the pumping in the convection 
zone and the depth at which pumping ceases (r 7 ) are un- 
known, we study the impact of varying these parameters in 
Section El 



2.3. Source term 

The source of poloidal field in this model is based on the 
observed emergence of bipolar sunspot groups. We use the 
procedure of CJSS10 to convert the RGO sunspot group 
area observations into a change in the radial magnetic flux 
at the surface: the RGO sunspot area record is used to ob- 
tain sunspot group areas, longitudes, and latitudes. Tilt an- 
gles from the Mount Wilson and Kodaikanal observa tories 
(|Howard et~aDll984l: lHowardlll99ll: ISivaraman et all 



1993) 



2010) 



2.2. Flow field and turbulent difFusivity 

We use the same profile for th e differential rotation as was 
used in ICameron et all ([20121 ) and the references therein. 
We also use the same form for the meridional flow as given 
in that study, except that we have reduced the parameter 
vq that determines the speed of the meridional flow so that 
the maximum speed at the surface is now 11 m/s instead 
of 15 m/s to be consistent with CJSS10. 



and CJSS10, and then used in conjunction with the RGO 
data for individual spot groups to determine the longitudi- 
nal and latitudinal separation between the two polarities in 
the bipolar group. The areas are converted into fluxes using 
the parameters in CJSS10 (also see the references therein). 
For each sunspot group we calculate the corresponding 
change in the radial field at the surface, AB r (R Q ,9, </>, t). 

The RGO data covers the period from 1874 to 1976. 
The tilt angles are known from 1906 to 1987, and the over- 
lapped data (cycles 14 to 20) can be used to reconstruct 
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Fig. 1. Source term, S(r, 9, t) sin(6>) at r = R® for the poloidal flux based on the RGO sunspot record. The factor of sin 9 
is included to make the figure proportional to the amount of signed flux. The red and blue colors indicate positive and 
negative signs of A, respectively. Data after 1977 is from the SOON network. 



the poloidal source term for cycles 15 to 21. The RGO 
sunspot records were replaced by results from the SOON 
network after this date, which used a different definition of 
what to include, and therefore introd uces a number of cross- 
calibr ation issues. As pointed out in ICameron fc Schiisslerl 
(|2012h there are outstanding problems when the sunspot 
record is used in a non-linear way such as here (the sunspot 
areas affect the tilt, which is then multiplied by the area to 
get the poloidal source term) . We thus restrict our calcula- 
tion to the period from 1874 to 1981 (i.e., for the cycles for 
which the toroidal flux at the base of the convection zone 
is created from poloidal flux based on the RGO area data). 

The tilt angle data is only available from 1906 onwards, 
so it can be used for the sources leading to cycles 15 to 21. 
For cycles 12 and 13, we use the linear relationship between 
the sunsp ot maxima and the cycle- averaged tilt angles re- 
ported bv iDasi-Espuig et al.l (|2010l) . While we restrict our 
statistical analysis to cycles 15 to 21, simulating these early 
cycles makes us less sensitive to the initial condition. 

To calculate the change in the azimuthal component of 
the vector potential A, we then use 



S(r,6,t)sml 



R Q f(r) 



2tt 



271 



'o Jo 

x 5{t-t cm ) (9) 

where t cm is the time at which the sunspot group was ob- 
served to emerge, and 5 is the Dirac delta function. For the 
radial dependence of the source term we take 



f(r) 



1 +erf 



O.9i? ) 



(10) 



The resulting sources at the surface are shown in Figure [T] 
This procedure is conc eptually related to the double -ring 
algorithm developed bv iMunoz- Jaramillo et al.l (|2010h . 

2.4. Initial conditions 

Since we have neither direct measurements of the magnetic 
fields at the poles of the Sun in 1874 nor measurements of 
the toroidal field of the Sun at the base of the convection 
zone at any time, the appropriate initial condition for our 



model is not well known. The starting date for our simula- 
tions, May 1874, is however about three fifths of the way 
through cycle 11 (which ran from March 1867 to December 
1878). This is close to the phase of the cycle where the axial- 
dipole field is weak, and given that the minimum of the ope n 
flux was low during this period (jWang fc Sheelevl 12009). 
A = is a reasonable choice. The poloidal component of 
the initial condition is particularly important because the 
decay of the poloidal field is dominated by the rate at which 
it is removed from the surface. The process by which the 
removal takes place is diffusion across the equator, which 
is opposed at the surface by the poleward meridional flow. 
The resulting decay rate is about 4000 years (see CJSS10) 
in the absence of new sources, so the poloidal component 
of the initial condition affects the entire period studied in 
this paper. 

The choice for the initial toroidal field is less critical 
because its diffusion rate is determined by diffusion across 
the equator at the base of the convection zone. Here the 
meridional flow is towards the equator, so it pushes the 
flux into a thin boundary layer where the field can easily 
diffuse across the equator. The pumping plays an impor- 
tant role here in that it keeps the meridional circulation 
from moving the flux upwards and back into the bulk of 
the convection zone. For this study we assumed B = at 
the start of the simulation. This may not be adequate, but 
in view of the short decay time of the toroidal field it is 
sufficient to demonstrate that the Babcock-Leighton model 
has polar fields at the preceding minima and toroidal fields 
at maxima that are strongly correlated with the maxima of 
solar activity. 

In summary, both the initial toroidal and the initial 
poloidal field arc zero in our simulations. We use a poloidal 
source term constructed from the observed sunspot record 
rather than the self-consistent dynamo problem. 



3. Results 

The time evolution of the simulated magnetic field through 
cycle 19 is shown in Figure [2] as an example. We see 
that throughout the cycle both the poloidal and toroidal 
fluxes have a simple form with no 'conveyor belt' carry - 
ing flux from multiple cycles (c.f. iDikpati fc Gilmani r2006'l 
oper ating within the convec tion zone, which is consistent 
with lKarak fc Nandvl (|2012fl . Although we are in the 'low- 
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medium diffusivity' regime (c.f. lYeates eta l. 2008), the 
transport by convective pumping of poloidal flux from the 
surface to the tachocline occurs over the same time span 
as the flux is transported to the poles. The winding up 
of the field to produce toroidal flux due to the differen- 
tial rotation takes place as the field is descending and after 
it reaches the base of the convection zone. The transport 
by downward pumping near the equator is opposed by the 
upward radial flow of the meridional circulation, together 
with the equatorward meridional flow, leading to a slower 
downward transport near the equator and a delayed rever- 
sal of the toroidal flux near the equator with respect to 
mid-latitudes. 

The radial component of the field on the surface as a 
function of time for the entire period analyzed is shown in 
Figure |31 The surface evolution of the field is similar to 
the one given by the surface flux transport model. Since 
with pumping the SFT and FTD models are consistent, 
and as we are using the same sources as in CJSS10, the 
correlation between the polar fields (here the average field 
from ±70° latitude poleward) and the strength of the next 
cycle (the maxima of the yearly averaged sunspot number) , 
based on cycles 15 to 21, is similar to the result in CJSS10 
(correlation coefficient r — 0.85). 

With our 2-D dynamo model, we can now also consider 
the evolution of the toroidal field at the base of the convec- 
tion zone, which could not be studied by the SFT model 
alone. Figure |4] shows the evolution of the toroidal flux at 
the base of the convection zone. The toroidal component 
of the magnetic field at r = 0.7 Rq is shown in the upper 
panel. In each hemisphere we see toroidal flux propagating 
equatorward at low latitudes. The reason for this propaga- 
tion is, in part, the delayed reversal of the toroidal fields 
there as discussed above. The eruption of this flux will re- 
sult in sunspot groups with an emergence latitude that mi- 
grates equatorward - that is, with properties similar to the 
observed butterfly diagram. There is also a high-latitude 
branch of toroidal magnetic field. We conjecture that this 
branch does not lead to sunspot emergence because, at lat- 
itudes above 75°, flux tubes at the bottom of the overshoot 
layer with field strengths of 6 x 10 4 G or more are stable with 
regard to the Parker inst ability (|Ferriz-Mas fe Schusslerl 
119951: ICalieariet all Il995h . This flux thus remains hid- 
den at the base of the c onvection zone. This explan ation 
contrasts with tho s e of [Nandv &; Choudhuril ([2002) and 
iGuerrero fe Munoa (|2004h , who show how a meridional flow 
that penetrates below the tachocline substantially weakens 
the high-latitude branch and pushes the toroidal flux below 
the tachocline. Flux emergence then can take place at low 
latitudes where the meridional circulation causes the flux 
to re-enter the convection zone. 

The evolution of the unsigned toroidal field is shown in 
the lower panel of Figure 4, taking proper account of the 
overlapping cycles at the base of the convection zone, at r — 
0.7Rq and averaged over latitudes -45° < A < 45°. The 
maxima of the modeled activity levels are clearly related to 
the observed amplitude of the cycle, with the correlation 
coefficient calculated using cycles 15 to 21 being r = 0.93. 

4. Effect of the pumping and diffusivity parameters 
in the bulk of the convection zone 

The dynamics in the bulk of the convection zone and near 
the tachocline are much less well constrained by observa- 



tions than are the near-surface dynamics. The strength of 
the pumping in the bulk of the convection zone (7 CZ ) and 
the turbulent diffusivity in the convection zone (j] C z) are 
therefore free parameters of the model. The evolution of 
the surface field is largely independent of these choices, but 
the flux at the base of the convection zone is expected to 
be affected. To investigate this we performed simulations 
for which we varied these two parameters, one at a time. 
The effect on the correlation between the toroidal field at 
the base of the convection zone and the level of activity is 
shown in Figure [5j The results show that the good corre- 
lation between the toroidal field at the base of the convec- 
tion zone and the activity maxima is robust for convective 
pumping of 2 to 10 m/s and for diffusivities in the bulk of 
the convection zone of 20 km 2 s _1 or less. 



5. Conclusions 

Our aim was to investigate whether the Babcock-Leighton 
model, with sources and parameters based on observations, 
is consistent with the correlations implied by the observa- 
tions. We used a flux transport dynamo model with pump- 
ing and a poloidal source term based on the historical record 
of sunspot groups' areas, locations, and tilt angles to study 
solar cycles 15 to 21, after starting the simulations during 
cycle 11. 

Following ICameron et al.l (|2012f ) , the inclusion of suffi- 
cient radial pumping of the magnetic field means that the 
surface field from the flux transport dynamo model, in par- 
ticular the polar fields, is similar to the surface flux trans- 
port results in CJSS10. In the present work, we have also 
been able to study the evolution of the subsurface field, 
including the toroidal field at the base of the convection 
zone. 

We found that the evolution corresponds to a particu- 
larly simple Babcock-Leighton type scenario: in each cycle 
the poloidal field from the surface is advected and diffused 
downwards so that by the end of a cycle, it completely 
replaces the previous poloidal flux. The toroidal field then 
results from the winding up of this poloida l field. The mem- 
ory o f the toroidal field is only one cycle (|Karak fc Nandvl 
2012). The decay time of the poloidal field is, in con- 
trast about 4000 years. The equatorward propagation of 
the activity belt is reproduced. We find that the simulated 
toroidal flux at the base of the convection zone is highly 
correlated with the observed sunspot area observed during 
the same cycle, consistent with its being the source of the 
emerging sunspots. The simulated minimum of the polar 
fields between cycles n — 1 and n are highly correlated with 
the maximum activity of cycle n. The Babcock-Leighton 
model thus appears to be consistent with and able to ex- 
plain the observationally inferred correlations. 
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Fig. 2. Simulated variation of the magnetic field over cycle 19. The red and blue colors correspond to negative and 
positive toroidal fields, respectively. The solid (dashed) lines represent anti-clockwise (clockwise) orientated held lines of 
the poloidal field. The snapshots (A)-(F) cover cycle 19, corresponding to t = 1956, 1958, 1960, 1962, 1964, and 1966, 
respectively. 
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Fig. 3. Radial component of the surface magnetic field for the reference case. The upper panel shows a time-latitude 
diagram. The lower panel gives the polar fields, defined as the average radial field strength beyond ±70° latitude (black 
curve) and the sunspot numbers (red curve). The vertical line indicates the time until which RGO sunspot group area 
data are available. The correlation coefficient of the maxima of the polar fields and the subsequent maxima of sunspot 
number for cycles 15 to 21 is 0.85. 
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Fig. 4. Toroidal field at the base of the convection zone (r = 0.7i?©). The upper panel shows the evolution as a function 
of latitude and time. In the lower panel the average unsigned toroidal flux between ±45° latitude corresponding to odd 
and even cycles is shown using dashed and solid curves, respectively. The correlation coefficient between the maxima of 
the toroidal field at the base of the convection zone and the maxima of the observed sunspot number for cycles 15-21 is 
0.93. 
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Fig. 5. Correlation coefficient between the maxima of the unsigned toroidal field corresponding to each cycle at the base 
of convection zone (r = 0.7R@) and the observed maxima of the sunspot number for simulations with different values of 
the pumping (left panel), and the turbulent diffusivity (right panel) in the convection zone. The asterisk represents the 
actual numerical experiments which were performed. The square indicates the result for the reference case. 
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